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Abstract. Quantum chromodynamics has a rather complicated phase structure. 
The finite temperature, chiral phase structure depends on the number of flavours 

Qand to a large extent on the particular values of the fermion masses. For two massless 
' flavours there is a true second order transition. It has been argued that this tran- 

\^ ' sition belongs to the universality class of the three-dimensional 0(4) spin model. 

The arguments have been questioned recently, and the transition was claimed to 
be mean field behaved. 

In this lecture we discuss this issue at the example of the three-dimensional, 
parity symmetric Gross-Neveu model at finite temperature, with a large number 
N of fermions. At zero temperature there is a phase where parity is spontaneously 
broken. At finite temperature, this model has a parity restoring second order tran- 
sition. It reveals considerable similarity to the QCD chiral phase transition. There 
are related questions here concerning the universality class. We solve this problem 
r~^ , essentially by means of the following methods: Large A'^ expansion, dimensional re- 

^\ • duction in the framework of quantum field theory, and high order convergent series 

expansions about disordered lattice systems. 

p i' 1 Introduction - QCD and the chiral phase transition 

D ' 

Quantum field theories provide the appropriate framework for the under- 
standing and the quantitative description of high energy physics. In particu- 
1^ I lar, quantum chromodynamics - QCD - is supposed to be the theory of strong 

'Oj . interactions, the physics of quarks, mesons and hadrons. As such a theory it 

C^ ' has to describe a large variety of high energy phenomena. For instance this 

concerns the confinement of quarks, the observed mesonic and hadronic par- 
ticle mass spectrum, topological excitations. This theory evolves an enormous 
complexity. Computational techniques to be developed and used are rather 
complicated, both numerical and analytical ones. 

There are a number of long-term projects running, with the aim of a 
quantitative understanding of particular properties of QCD from first prin- 
ciples. For example, one of these projects is the so-called alpha-collaboration 
(Liischer, Sint, Sommer, Weisz and Wolff (1996)). One of the aims is to deter- 
mine the flow of the (appropriately defined) renormalized coupling constant 

* Supported by a Heisenberg Fellowship 



2 Thomas Reisz 

Qr^L) well above the perturbative short distance scale, that is on large spa- 
tial scale L. A remarkable strategy of this collaboration is the combination 
of the appropriate analytical and numerical (Monte-Carlo) techniques. 

Other long-term projects concern the nature of the chiral phase transition 
in finite temperature QCD and the properties of the high temperature phase. 
The high temperature phase of QCD may be realized in heavy ion collisions, 
and it was certainly realized in the early universe. There are a large number 
of interesting physical questions related to this phase. For instance, screen- 
ing of heavy quarks in the QCD plasma phase has been studied in detail 
and screening masses have been extracted, by combining analytical methods 
such as finite volume perturbation theory and dimensional reduction in the 
framework of quantum field theory, together with various Monte-Carlo meth- 
ods. Other, yet open questions concern the (gauge-invariant) definition of a 
magnetic mass and its determination, the ability to disentangle infrared and 
ultraviolet singularities and their removal in the thermodynamic limit, that 
is to find appropriate resummation techniques. 

In the limit of vanishing current quark masses, massless QCD has a chiral 
symmetry. This global symmetry is spontaneously broken at zero tempera- 
ture, and it becomes restored in the high temperature phase. Whether this 
symmetry is restored by a 1st or 2nd order phase transtion or just by a 
crossover phenomenon depends on the number of fermionic flavours. In ad- 
dition, to a large extent it depends on the particular values of the fermion 
masses. 

A detailed presentation of the QCD phase transitions according to our 
current knowledge can be found in the review of Meyer-Ortmanns (1996). 
Here we only give a very short outline of the state of the art. In Fig. 1 we 
show a qualitative plot of the supposed phase diagram for three flavours u, d, s 
in dependence on their masses. The u- and d-quark masses are of comparable 
size and arc put to the same value. For the pure SU(3) Yang-Mills theory 
(right upper corner) the transition is the deconflnement transition, which is 
believed to be of 1st order. For the purely massless case (left lower corner) 
the transition refers to the chiral transition, which is supposed to be of first 
order in the case of 3 flavours. With increasing masses the chiral transition 
becomes weaker and finally turns over into a mere crossover. For 2 massless 
flavours, the chiral transition is of first order below a critical value m* of the 
strange quark mass (corresponding to a tricritical point) and of second order 
above to*. 

The physical fermion masses are subject to a large uncertainty in relation 
to the phase boundaries. Even if there is no "true" phase transition in QCD, 
some correlation lengths might become large, generating large domains of 
chiral condensates. 

In nature, two fermion flavours are realized as relatively light states. For 
the case of two flavours, both being massless, the chiral phase transition 
is supposed to be described by a critical point. What is the nature of this 
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Fig. 1. A qualitative plot of the mass dependence of the QCD phase transition for 
three flavours. The dot identifies the tricitical strange quark mass, above which the 
transition is of 2nd order. 



transition? It was argued by Pisarski and Wilczek that its universality class 
is that of the 0(4) spin model in 3 dimensions. We will discuss this point of 
view in the next section. On the contrary, Kocic and Kogut argued that the 
transition is mean field behaved. So far, no definite answer has been given, 
neither by Monte Carlo simulations nor by analytical investigations. 

In this lecture we discuss the issue of the finite temperature phase tran- 
sition at the example of the 3-dimensional Gross-Neveu model. This model 
reveals considerable similarity to QCD. There are quite related questions con- 
cerning the nature of the transition, including as candidates both mean field 
and non-trivial spin model behaviour. However, in this case a definite answer 
can be given, mainly by analytical means. 



1.1 Two-flavour QCD 

Before we turn over to the Gross-Neveu model, let us discuss in some more 
detail the effective scalar 0(4) model for the QCD chiral phase transition. 
The order parameter is given by the chiral condensate < ip{x)ip{x) >. It 
is argued that the transition belongs to the universality class of the scalar 
0(4) model in 3 dimensions (Pisarski and Wilczek (1984), Wilczek (1992), 
Rajagopal and Wilczek (1993)). An essential ingredient of this argument is 
the observation that under quite general conditions, for a field theory at high 
temperature fermionic degrees of freedom act only as composite states in the 
infrared. This fact is a consequence of the so-called dimensional reduction in 
quantum field theory and will be discussed in detail in the following sections. 
Let us outline the way of reasoning towards the 0(4) scenario. The idea 
is to build an effective model of the "localized" order parameter field, the 
condensate, and apply a Landau- Ginzburg type of argument. First of all, 
the order parameter field is a complex 2 x 2-matrix field in three dimensions 
X gR3, 
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< ^mix)^^Lj{x) >^ My(x) e M(2 X 2,C), (1) 

where R and L denotes right- and left-handed (relativistic) fermions, respec- 
tively. Close to the critical point the degrees of freedom that develop long- 
range fluctuations are assumed to be associated with the order parameter M , 
which is small in magnitude. Hence it is plausible to describe the effective 
model by the action 

Seff{M)= f {tTd^M{x)''d^M{x) + m^trM''{x)M{x) (2) 

J X 

+31 tr (Mt(a;)M(a;))2 + g2itT M\x)M{x)f). 
The partition function becomes 

Z= fvM -ex^pi-SeffiM)), (3) 

where VM — Yix'^-^i^)- ^^ have implicitly assumed an ultraviolet cutoff, 
such as the Pauli-Villar cutoff or a lattice. At this stage then we can exploit 
the methods of renormalizable quantum field theory, such as the renormal- 
ization group and renormalized perturbation theory, and strong statements 
obtained by those methods. On the one hand, a posteriori this justifies a 
large freedom in the choice of the ultraviolet cutoff for the above (super-) 
renormalizable model. That is, a particular choice does not have influence on 
the low-momentum scale. On the other hand, the critical properties of the 
model can be studied. 

Yet the model has too large a symmetry to describe a chiral phase tran- 
sition. It is invariant under global U{2)l x U(2)ji transformations 

M ^ UMV'' , U,VeU{2). (4) 

Finally we want to reduce the symmetry to the smaller, physical symmetry 
group SU{2)l X SU{2)r x U{1), (to be broken down to vector SU{2) x U{1) 
spontaneously at low temperature,) where U{1) describes the vector-baryon 
symmetry. This restricted symmetry has to be implemented by imposing 
appropriate constraints on M . Whereas for larger number of flavours those 
constraints are to be non- linear in M, leaving the framework of renormalizable 
quantum field theory, for two flavours we can implement such constraints by 
restricting M to be a real multiple of a unitary matrix of unit determinant, 
M ^ c ■ SU{2) , c e R. Hence M can be parametrized linearly as 

M = 7rol2 + 7rcr, (ttq, tt) e R'', (5) 

where cr = {cri,a2,<J3) denotes the Pauli-matrices. It is now obvious that in 
terms of the 7r-fields we get a linear 0(4) model. This symmetry is broken at 
low temperature to 0(3). 

If the universality arguments made so far are true, the QCD chiral phase 
transition for two massless flavours will belong to the universality class of 
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the 0(4) spin model in three dimensions. Of course, there are assumptions 
made in the course towards this conclusion. The most stringent ones are the 
following. 

— The chiral phase transition is of 2nd order. 

— The free energy density as a function of the order parameter field M is 
sufficiently smooth, allowing for a polynomial representation at small M . 
But we know this is problematic below the upper critical dimension, the 
dimension of strict renormalizability, which is 4 for the above model. 

— In the language of statistical mechanics, the quantum critical system is 
equivalent to the corresponding classical critical system. In particular, 
the phase transition temperature is so high that the thermal fluctuations 
dominate the quantum fluctuations even at the transition. In the language 
of quantum field theory, this amounts to say that dimensional reduction 
is complete at the critical point, at least concerning the values of the 
critical exponents. 

Recently Kocic and Kogut have criticised the above (a-model) approach 
because the fermionic degrees of freedom are represented only via composite 
fields M (Kocic and Kogut (1994)). Actually, as we will see in the course of 
this lecture, this is a minor point of criticism. But even if fermions act only as 
composite states, there are several candidates for the universality class of the 
QCD chiral phase transition. Both the above argumentations mainly refer to 
conventional wisdom. In order to obtain a definite answer it is unavoidable 
to give strong reasonings, including explicit computations. 

1.2 Why studying the Gross-Neveu model? 

In this lecture the problem of the chiral phase transition will be discussed at 
the example of a much simpler model than QCD, which nevertheless bears 
the essential properties such as " chirality" . Namely, I would like to discuss the 
3-dimensional Gross-Neveu model at finite temperature, with a large number 
of fermionic flavours. It is a model with a four-fermion interaction, subject 
to a global Z2 symmetry (which is parity). This symmetry is known to be 
broken at zero temperature, and it is restored at high temperature by a 2nd 
order transition. We then ask for the universality class of this transition. 

In quite analogy to QCD there are two opposite claims of the universal- 
ity class. Adjusting the arguments of Pisarski and Wilczek, the transition is 
described by a Z2 symmetric scalar model in 2 dimensions that should be- 
long to the Ising class. On the contrary, Kocic and Kogut claimed that the 
transition is that of mean field theory. 

Beyond its similarity to QCD, the advantage of studying this model is 
that to a large extent we have a closed analytic control of it. This includes 
the finite temperature transition. As we will see, the above discrepancy can 
be resolved by analytical means. Also, in the course of the following sections 
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it will become clear that the definite answer can be given only by explicit 
computations. 

The answer will be given essentially by applying the following two tech- 
niques. 

— Dimensional reduction in quantum field theory. Applied to the 3-dimen- 
sional, finite temperature Gross-Neveu model, it states that the infrared 
behaviour of the latter is described by a 2-dimensional scalar field the- 
ory at zero temperature. As we will see, this effective model is local and 
renormalizable. It should be emphasized that dimensional reduction pro- 
vides an explicit computational device to compute the couplings of the 
effective model. This is important because a 2-dimensional scalar model 
does not necessarily belong to the Ising universality class. 

— Convergent series expansions of high orders, applied to the free energy 
and connected correlations. In this way the critical region can be inves- 
tigated, the order of the transition and eventually the universality class 
can be determined. 

In outline, this lecture is organized as follows. Throughout we will use 
the path integral representation for quantum field theories. In Sect. 2 we de- 
scribe the 3-dimensional Gross-Neveu model. The notions familiar from finite 
temperature quantum field theory are introduced. The large A^ representa- 
tion and expansion of the model are discussed. The N ~ oo limit can be 
solved in closed form, including the phase structure. In Sect. 3 we outline 
the technique of dimensional reduction in the framework of renormalizable 
quantum field theory. In the recent past this technique became a very well 
developed, powerful computational machinery both for bosonic and fermionic 
field theories at finite temperature. It is applied to the Gross-Neveu model. 
Then we are prepared to study in Sect. 4 the phase structure of the model 
for finite number A^ of flavours. A first insight is obtained by the discussion 
of the strong coupling limit. In a subsection we shortly present convergent 
series expansion techniques, in particular linked cluster expansions, and the 
ideas behind their generation. They are then applied to determine the criti- 
cal exponents, hence the universality class of the parity restoring transition. 
Summary and outlook for the case of QCD is given in the last Sect. 5. 

2 The 3-dimensional Gross-Neveu model at zero and 
finite temperature 

2.1 The model 

For a detailed introduction to field theories at finite temperature, we refer to 
the books of Kapusta (1989) and of Rothe (1997). We first of all define the 
model under consideration. In 3 dimensions at finite temperature T (hence- 
forth Sdr in short), a field theory lives in a volume of the toroidal shape 
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z = (zo,zi,Z2)eee(zo,z) eR/TxR2, (6) 

cf. Fig. 2. The extension is set to infinity in both spatial directions, and the 
length of the torus component, the "temperature-direction", is fixed by the 
inverse temperature. 



R^ 




Fig. 2. The world of the Sdr Gross-Neveu model. The fermionic degrees of freedom 
are subject to anti-periodic boundary conditions in the temperature direction. 



We consider a model of N Dirac spinor fields, that is of N fcrmion species 
or flavours, with values in a Grassmann (exterior) algebra, 

-tp{z) = {i'aiiz))\a=l,2;t=l.,...,N, 

^{Z) = {i^o.i{z))\o.= 1.2■,^=l,■■■.N■ (7) 

Due to the trace operation of the partition function, Grassmann- valued fields 
are subject to anti-periodic boundary conditions in the temperature direction, 

i/.(zo±T-i,z)--^(zo,z), (8) 

and similarly for ^/) (cf. e.g. Liischer (1977)). As a first consequence of the 
anti-periodicity, the fermionic degrees of freedom always have an implicit 
infrared cutoff involved, which is of the order of the temperature T . This is 
easily seen by applying a Fourier transform to momentum space, 

V^(z)= E e^'"'" ^«(z). (9) 

fco=7rT(2n+l),neZ 

From this we infer that |fco| > ttT, which provides a lower bound on the total 
momentum. 

The action is given by 

Sgni^A) = j(^{z)D^(z)) - ^ ^ {^{zmz)f , (10) 

where 

1 rp-l 

r r „ /■ 2 ^ 

dzQ, 



JR2 




2 


d 



D = l^l^—, (11) 

i=0 
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with 7o , 7i , 72 denoting the 3 basis elements of a Dirac representation. We 
normalize them according to 

l^l] + irn = 2<5»j ; i, J - 0, 1, 2. (12) 

Here and in the following for spinor fields tp{z)^ 77(2) we use the notations 

tpiz)r]{z) = Y^ Y^ Vai(2)?7m(z). (13) 

Q = l,2 i=l,...,N 

The partition function (generating functional of full correlation functions) is 
then given by 



Z{r^,rj)^ J Vi,V4> exp(-5<,„(^,^) + y (7/^(z)r;(z)+^(z)V(z))), (14) 

where 

jvi^V^lj = II d4,a,,{z)(hP^,iz). (15) 

7] and 77 are external sources introduced in order to derive correlation functions 
from Z{r],rf) by differentiation. 

So far, the model is only formally defined by the action (10) and the par- 
tition function (14). In order to define it rigorously as a quantum field theory 
from the outset, we have to introduce an ultraviolet cutoff A. Among others, 
possible choices are a lattice, the Pauli-Villar cutoff, dimensional regulariza- 
tion, or just a simple momentum cutoff. The model becomes well defined for 
p < A, where p represents momentum. An essential point then is to prove 
the existence of this model in the large cutoff limit yl ^ 00, and in this limit 
a sequence of axioms have to be satisfied (Osterwalder and Schrader (1975), 
Zinoviev (1995)). This is the issue of renormalization in the framework of 
quantum field theory. In many cases this is achieved by adjusting a finite 
number of bare parameters of the action, fixing corresponding renormalized 
coupling constants on a finite length scale. Renormalizability implies that the 
infiuence of the cutoff A on the low-lying momentum scale p is suppressed by 
some power of p/ A. In particular, all the correlation functions of the model 
stay well-defined in the limit on finite scales. 

One might argue that in the framework of statistical mechanics, one is not 
necessarily interested in the removal of the ultraviolet cutoff. For instance, 
the lattice spacing a ^ 1/ A is a finite physical length in sohd state physics 
(except for the lattice spacing in the temperature direction, which still is to 
be sent to zero). However, concerning the critical regions of a model where 
the correlation length diverges in units of the lattice spacing, the motivation 
for renormalizability is quite similar as in quantum field theory. Instead of 
the cutoff being sent to infinity, the range of interest is where the momentum 
becomes arbitrarily small. Collective critical behaviour in this region to a 
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large extent becomes independent of the cutoff scale, i.e. of the bare param- 
eters. This identifies a universality class, that is universal critical behaviour. 
This would not be true for non-renormalizable models. 

In the following we regularize the Gross-Neveu model by choosing a sim- 
ple momentum cutoff A. This is done mainly for computational simplicity. 
Equally well we could have chosen a hybercubic lattice, with some care con- 
cerning the choice of lattice fermions. 

There is some detailed knowledge of this model at zero temperature 
T — 0. First of all, applying power counting of ultraviolet divergencies, this 
3-dimensional model appears to be non-renormalizable in the loop expansion. 
That is, as an expansion in powers of the coupling constant A, we need more 
and more couterterms with increasing order of A in order to render the ultra- 
violet limit finite. Fortunately, for large number N of fermionic fiavours, the 
situation is different. Within the large N expansion, which will be explained in 
the very detail in the next section, the model becomes renormalizable, to ev- 
ery finite order in 1/A^ (Park, Rosenstein, Warr (1989a)). Actually, there is an 
even stronger statement. Namely, if N is sufficiently large, the 3-dimensional 
Gross-Neveu model can be shown to exist in the large cutoff limit (and in 
the infinite volume), that is, the correlation functions of the model stay fi- 
nite. The proof is done in the framework of constructive field theory, using 
convergent expansions and the renormalization group (Faria da Veiga (1991), 
de Calan, Faria da Veiga, Magnen, Seneor (1991)). Remarkably enough, this 
is a construction about a non-Gaussian fixed point, leading to a non-trivial 
(interacting) quantum field theory. 

Let us come back to the finite temperature. As already said in the in- 
troduction, the model (10) serves as a good example for studying a "chiral" 
phase transition because it reveals essential similarities to massless QCD. We 
seek for a symmetry that would be broken explicitly by a fermionic mass 
term of the form 

?Til/'(z)'0(z). (16) 

Actually, this will not be a true chiral symmetry. This would require the exis- 
tence of a 7- matrix, the 75 , which anticommutes with all generating elements 
(7o:7i7 72) of a faithful representation of the associated Cfifford algebra (cf. 
e.g. Chevalley (1997) or Porteous (1995)). It does not exist in three dimen- 
sions. Instead, the model is invariant under the global parity transformation 

ip{z) -^ ip{-z), 

^(z) ^ -7j{-z). (17) 

Both the action (10) and the measure (15) stay invariant under this trans- 
formation. The symmetry under (17) is spontaneously broken at T = 0, and 
it becomes restored at high temperature by a second order phase transition, 
cf. below. The parity condensate 
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- ^(^ / ^^^^ ^ (^(^) V'(^)) exp (-5,„ (?, V)) (18) 

serves as an order parameter of the transition. 

We remark in passing that, besides others, the model has a continuous 
global symmetry in flavour space, defined by 

N N 

tpaiiz) -^ Yl Uijlpajiz) , IPatiz) -^ Yl ^li^ajiz), (19) 

with U G U{N). This symmetry is neither explicitly broken by a mass term 
(16) nor does it get broken spontaneously. In particular it yields fermion 
number conservation. 

2.2 Many flavours. The large TV expansion 

As mentioned in the previous section, much insight is obtained by represent- 
ing the Gross-Neveu model in a form that is particularly well suited if the 
number of fermionic species is sufflciently large. The corresponding transfor- 
mation is based on the resummation of so-called planar graphs according to 
Fig. 3. 



+ 



Fig. 3. Resummation of planar graphs. 
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Here c counts the combinatorical factor and the trace in spinor space over 
the single fermion bubble. Such a chain then is represented as an interaction 
mediated by an auxiliary scalar field a. 

We can do this resummation in closed form in the following way. First 
the local, quartic fermion interaction part of the action can be written as the 
result of a gaussian integration, using 



e.,[^mm') - (j^)' 



1 / N 



da{z) exp ( -- ( — ) a{zY ± a{z) i^{z)^{z) ] . (21) 
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The action then becomes of the Yukawa-type and is biUnear in the fermionic 
fields. In turn, the integration over the Grassmann algebra valued fields can 
be done according to 



DV^V exp ( - / / V(2) Q{z, z')'ip{z') + / {xp{z)r){z) + rj{z)iP{z)) 

= det(-Q) exp n j^^ v{z)Q-\z, z')7j{z')) . (22) 

In our case, Q is flavour diagonal, i.e. Q = 1 (g) K, with 

K = D + a-l (23) 

acting on the product of spinor and configuration space. D is defined by (11). 
As the remnant of the anti-periodic boundary conditions (8) imposed on the 
fermionic fields, the momentum space representation of D has nonvanishing 
components only for the discrete thermal momenta fco ~ 7rr(2m + 1), m G Z, 
cf. (9). In this way we obtain the following representation for the partition 
function (14) up to an unimportant normalization constant. 

Z{ri,rj) = const / Va exp {—N Seff{cr)) 

■ exp (J^ J^^ rj{z){l ® K-^)iz, z')rj{z')) , (24) 



with 



fva = Yldcr{z). (25) 



The efl:ective action 5*6// (cr) is given by 

Seffi^) = ^ (^ - tr, (logX) (z,z)) , (26) 

with the trace taken over the spinor indices. A careful counting of the degrees 
of freedom (e.g. by imposing a lattice cutofi^) shows that without loss of 
generality we may impose periodic boundary conditions on the auxiliary field 
a{z) in the temperature direction (it is "bosonic"). This is the natural and 
common choice and implies that constant field configurations of a are to be 
identified as zero momentum configurations in the usual way. In turn the 
minima of the effective action for large N belong to the latter ones. 

It should be emphasized that up to now no approximation is involved. 
The partition function has only been rewritten in a form that is most conve- 
nient if the number TV of fermions becomes large. The only place where the 
number A^ enters in (24) is as the prefactor of the action. The explicit form 
(24) reveals the standard saddle point expansion as the convenient method 
to study the Gross-Neveu model for large A^. The parity condensate (18) 
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becomes the expectation value of —ti sK^^{z, z), and similar for the conden- 
sate correlations. It is straightforward to show to all orders of 1/iV that a 
vanishing condensate (18) is equivalent to a vanishing expectation value of 
the auxiliary field a. The auxiliary field <j{z) itself is parity odd. The parity 
symmetry of the originial model becomes invariancc under 

a{z) ^ -a{-z). (27) 

The large N expansion proceeds along the following steps. 

— Localization of the translation invariant minima of the effective action 
Seff- This yields the so-called gap equation 

^^^U)^,^0\ (28) 

— Minimization requires positivity of 



'I J I 



5a{z)5c 
Expansion about the solutions according to 



(29) 



<r{^) = M + ^. (30) 

Existence of stable solutions of the gap equation and their properties will 
be discussed in the next section. For further reference we write down the 
model as obtained by inserting the ansatz (30) into the effective action (26) 
and expand for large N . Towards this end we first introduce some notations 
used in the following. Fermionic and bosonic thermal momenta are defined 
as the sets 

T={^{2m+l) |meZ}, 

P 

B={^2m\m^Z]. (31) 

For the inverse temperature we write (3 — T^^ . For x, /3 > and n = 0, 1, 2, . . . 
let 



Jn{x,(i) = 7^ 51 



' d2q 1 



a ^ J (27r)2 [q^+xY' 
got./" 

The prime at the integral sign denotes that the sum and integral is confined 
to the region with 

g2=g2+q2 <^2^ (33) 
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For /3 < oo the J,i(x, /3) are always infrared finite for every x > 0, whereas 

J„(x,^)=/ ^{q^ + x)- (34) 

is infrared singular as x ^ 0+ for n > 2. Furthermore, the ultraviolet singu- 
larities of the J„ arc temperature independent, which means that 

lim (J„(x, /?) — J„(x, oo)) < oo, (35) 

for each x > 0. This is a property which is shared by a large class of 1-loop 
Feynman integrals. They can be rendered ultraviolet finite by imposing the 
appropriate zero temperature normalization conditions. 

Inserting (30) into the effective action (26), we obtain the following large 
A^ representation. Let 

Ko^D + fx. (36) 

The partition function then becomes 
Z{T^,rj) ^N jV(t> exp (- S{(f>) + j J ^^ ^(^)(1 ® K-^)i^, z'M^')) , (37) 



with 



A/" = const exp(-A^ f {AX^fi^{Ji{fi^, (3)f - trs\og{Ko){z,z))Y (38) 



K = K„ + -^, (39) 

N2 



and 

^2 



Si<P) = ^ (^ - N tr, (^log(l + -^cj^K,') - ^'^^o^') i^^^\ 

= J2SnW- (40) 



n>2 



Here, 



S2{4>)-l{^^ + ltr4cf>K^')\z,z)^ 



Sn{<t^) = j^^i - ltr,{<I^K-^T{z,z), n>3. (41) 



n 



Sn{4>) is of the order of 0(iV^^"/^). For later convenience we write the 
quadratic part of the action in its momentum space representation obtained 
by Fourier transform. With 

'^(^) = ^ E ^^'" ^(fc) (42) 



14 Thomas Reisz 

we obtain after some algebra 



koGB 






' W2 






/fc2 + 4^2 



1 1 



(43) 



{q + A:)2 + ^2 ^2 ^ ^2 

The wave function renorniaHzation constant Zr and the renormalized mass 
TO/{ (the inverse correlation length) are defined by the small momentum be- 
haviour of the (connected) 2-point function 

W^'^\k) ^<4,{k)4>{-k) > (44) 

according to 

W^^\ko^O,k) ^ ^(ml + k^ + 0{k*)) as k ^ 0. (45) 

Zr 

2.3 Phase structure in the AT = oo Hmit 

We discuss the phase structure of the Gross-Neveu model in the limit of 
infinite number of fermionic flavours first (Park, Rosenstein, Warr (1989b)). 
The gap equation (28) reads 

the solution of which yields the condensate fJ.{l3). We are mainly interested 
in the phase structure at finite temperature, but it is instructive to include 
the case P — oo oi zero temperature, /i = is always a solution of (46). For 
fixed (3, Ji(/x^,/3) is monotonically decreasing with ^^. Hence, for 1/(4A^) > 
Ji(0,/3), /x = is the only solution, and it is straightforward to show that 
it is a minimum of the action. On the other hand, for 1/(4A^) < Ji(0, /?) 
there are two non-trivial stable solution a{z) = ±/i(/3) of (28) with /x(/3) > 
satisfying (46). 

Now suppose that the bare coupling constant A is so large that the parity- 
broken phase is realized at zero temperature. We define a zero temperature 
mass scale M by 

Ji(M2,^)ee^< Ji(0,^). (47) 

For fixed /i^, as a function of f3, Ji(/i^,/3) is monotonically increasing in (3. 
Decreasing (3 down from (3 = oo, there is a unique critical value (3c where the 
broken phase becomes unstable, cf. Fig. 4. This (3c is given by 
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Ji(0,cx)) 
Ji(M2,oo) 



MO, (3) 



high T symmetric (/i — 0) 




Tc 



T = p- 



Fig. 4. Phase structure of the Sdr-dimensional Gross-Neveu model at A'^ = oo. For 
given mass scale M, the symmetric high temperature phase becomes unstable at 
Tc = M/{2 In 2), and the model undergoes a second order phase transition to the 
broken phase. 



Ji(m2,w) = Ji(0,/3,) 



/3c = 



2 In 2 
M ' 



(48) 



In order to show that the symmetry gets actually broken spontaneously 
for P > Pc, before the infrared cutoff is removed (the volume is sent to 
infinity) we introduce an explicit parity breaking term (16) into the original 
fermionic action (10). Without loss we take tti > 0. Then we remove the 
infrared cutoff and finally let m -^ 0+. The only modification the fermionic 
mass term implies for the derivation of the large N representation of the 
partition function as done in the last section is the replacement of (23) by 



K 



D + im + (T) -l. 



(49) 



The gap equation (46) is replaced by 



(fJ'-m) — -M(M^/3) = 0. 



(50) 



For small m and (5 > (ic there are two stable solutions of (50) given by 
/z(m, (3) = ±n{P) + S, where /i(/3) > is the "unperturbed" solution of (46) 
and 

'^8AW)v'(,(/3)^/3)™ + ^^™^■ ^''^ 

Inserting this into (38), we see that the negative solution of the gap equation 
gets exponentially suppressed in the thermodynamic limit for positive to. We 
are left with a unique absolut minimum of the action, breaking the parity 
symmetry. 

The phase transition is of second order. The correlation length m^ di- 
verges at f3c- We define critical exponents in the standard way, in particular 
with X2 as the 2-point susceptibility. 



X2 



1^(2) (fc 



0)-^. 


--{Pc- 


-/3)-^ 


TO_R 2 


--{(3c- 


-/3)^ 


Zr- 


^{Pc- 


-pr, 



as (3 ^ Pc 



(52) 
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and similarly primed exponents v' ,^' by approaching the phase transition 
from the broken phase. To leading order in N~^ the 2-point function is given 

by 

W?(2)(fcQ ^ o,k) = 2 (Ji(M2,oo) - Ji(0,/3)) 

+ k2 2(33(0, P) + 0(k4) (as k ^ 0) (53) 



(/3,-/3)+ k^-L^ + 0(k^) 



7r/?/3c ' ' 167r 

in the symmetric phase, and 

W?(2)(fco = 0,k)-V^2(A*',/?) 

+ k^ 2 (J2(m',/3) - 2^?M^l\p) + A^?Q^{^l\p) - g3(A^',/3)) 

+ ©(k"*) (as k -> 0) (54) 

^ ^ (/3 - /3c) + k^l^ + 0(k4) (as /3 ^ /?,-) 

in the broken phase. We thus obtain ^f = ^' = \, v = v' = 1/2 and 77 = 0. 
Other critical exponents are determined in the straightforward way, with the 
result that a = a' = 0, (5 = 3 and /3 = 1/2. 

2.4 What happens for finite AT? 

So far we know the phase structure of the ZcIt Gross-Ncveu model in the 
limit of infinite number of flavours. As next, we are interested in the physics 
of the model if there is a large but still finite number of fermionic species. 
The parity symmetry still remains broken at small temperature and will be 
restored at sufRciently high temperature. Concerning the universality class 
of the proposed second order transition two opposite claims have been made. 

On the one hand, one might follow the resoning proposed by Pisarski and 
Wilczek. The length of the temperature torus is given by the inverse temper- 
ature /3. In the high temperature regime this length is rather small, hence the 
geometry is essentially two-dimensional on length scales large compared to 13. 
Furthermore, due to the anti-periodic boundary conditions to be imposed on 
the fermionic degrees of freedom, the latter are subject to an infrared cutoff, 
cf. (9). Hence, in the far infrared, only the purely scalar modes survive. It 
might be natural to assume that this provides a model that belongs to the 
Ising universality class in two dimensions. This would imply the following 
values for the critical exponents: i^ = 1, 7 = 1.75, rj = 0.25. This is to be 
contrasted to the case of A^ = 00. 

On the contrary, the reasoning proposed by Kocic and Kogut roughly goes 
as follows. The model is subject to a Z2 symmetry. This symmetry does not 
depend on the number of fermionic flavours. For N = 00 the model is gaussian 
(y = 1/2, 7 — 1), hence it might be so for finite N. A basic assumption made 
in this way of reasoning is that it is the symmetry pattern which determines 
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the universality class of a phase transition. But we know that this can be 
completely misleading. For instance, let us consider a classical spin model 
in 3 dimensions as described by an action of the type (jP + A^**. It has a 
global symmetry (/) ^ —(f> that does not depend on a particular value of the 
quartic coupling A. For A = the model is purely Gaussian. But this is in 
sharp contrast to the critical behaviour at finite A. Similarly, for the Gross- 
Neveu model, the large N behaviour is not predicitve for finite N, at least 
not without closer inspection. 

Which scenario is the right one? The reduction argument as applied above 
is purely geometrical. Does it generalize to field theory, and if so, how does 
the effective 2 dimensional model look like? Renormalization group studies 
reveal an infinite number of fixed points of those models, likely to be related 
to various different universality classes. It is not just the symmetry pattern 
alone that determines the critical behaviour. To decide the nature of the 
transition requires more sophisticated methods known from quantum field 
theory. 

Below we shall combine two techniques to resolve the puzzle of the Gross- 
Neveu model. The first one is dimensional reduction in the framework of 
quantum field theories. This has become a well established machinery in the 
recent past. It provides a computational device to derive effective models for 
the quantitative description of the infrared behaviour at high temperature. 
The second method concerns convergent high temperature series. Under very 
general conditions, the high order behaviour of those series allows for the 
study of the critical region, including the measurement of critical exponents. 



3 Dimensional reduction in quantum field theory 

Field theories for realistic systems and statistical models with a large number 
of degrees of freedom evolve considerable complexity. From the point of view 
of high energy physics, good examples are quantum field theories such as 
quantum chromodynamics, which is supposed to be the microscopic theory of 
the strong interaction, and the electroweak standard model. The investigation 
of the large scale properties becomes very demandingful. In many instances 
we are forced to use methods that imply appropriate simplifications. On the 
one hand we will learn for the more complicated situation. On the other hand 
it is desirable to have a systematic approach, such as convergent or asymptotic 
expansions. This allows for a quantitative investigation and statements to be 
made for the full theory. 

For field theories at finite temperature, one of those systematic techniques 
is dimensional reduction. Originally it served as a more or less qualitative rea- 
soning for the description of the infrared behaviour at high temperature. It 
is easy to understand and works well for classical field theories. The general- 
ization to quantum field theories, however, is neither simple nor obvious. The 
first attempts were based on a naive generalization, by just assuming that 
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decoupling properties of heavy masses that are well known for zero temper- 
ature generalize to finite temperature. It took some time until it was made 
definite that this is not the case (Landsman (1989)). 

A systematic way in which dimensional reduction is realized in renor- 
malizable quantized field theories was given by the author (Reisz (1992)). 
Meanwhile it has become a very powerful technique. In particular it has been 
used to study the infrared properties of the QCD high temperature phase, 
so as to extract screening masses of the quark-gluon plasma (Karkkainen, 
Lacock, Petersson and Reisz (1993)). For a summary and further references, 
cf. Lacock and Reisz (1993). Later, the success of this well developed tech- 
nique and the experiences obtained inspired its application to the electroweak 
standard model of the early universe (Kajantie, Laine, Rummukainen and 
Shaposhnikov (1996)). In the following we will shortly outline the ideas of 
dimensional reduction and then will apply it to the Gross-Nevcu model. 

3.1 A short survey of dimensional reduction 

We consider a local, rcnormalizable quantum field theory in thermal equi- 
librium by the interaction with a heat bath, in £) > 3 dimensions. In the 
functional integral formalism such a field theory lives in a D-dimensional vol- 
ume with one dimension being compactified to a torus of length given by the 
inverse temperature (Fig. 2). There are boundary conditions imposed on the 
physical fields in this direction due to the trace operation of the partition 
function. Bosonic degrees of freedom are subject to periodic boundary con- 
ditions, fermionic degrees of freedom to anti-periodic boundary conditions. 

The idea behind dimensional reduction is based on the geometrical picture 
that at high temperature T the length of the temperature torus given by T~^ 
becomes small. So as the temperature increases, the D-dimensional system 
becomes more and more D — 1-dimensional on spatial length scales L large 
to T^^. In the following we discuss the implications for a rcnormalizable 
quantum field theory. 

Consider, for simplicity, a scalar model with field content 4'{z), described 
by the action 

SW^l dzoJd''-'z(^q^{z){-AD^){z) + V{cj,{z))Y (55) 

where Ad = X]i=o ^^/(^^?) denotes the D-dimensional Laplacian with peri- 
odic boundary conditions in zq, and V the self-interaction of the field. Period- 
icity implies that the corresponding momentum components in this direction 
are discrete. We apply a so-called Matsubara decomposition, which is nothing 
but a Fourier transform with respect to this momentum component, 

0(z) = tV2 ^ e^feo.o ^^(^) 

fco=2n7rT,neZ 

= Ti/2 0,,(z) + 0„,(z), (56) 
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where the static field (f>st denotes the n = component and the non-static 
field (f)„s sums all components with n ^ 0. The kinetic part of the action 
becomes 

^ j d""-^^ \ ?o(z) (-Ad-i4>o) (z) (57) 

+ Y. I d""-^^ \ ?„(z) (^(-Ad-iK) (z) + (2^nr)2^„(z)) . 

Apparently the field components with n ^ acquire a thermal mass pro- 
portional to T. This provides an infrared cutoff to these non-static modes 
mediated by the interaction with the heat bath. The only modes which do 
not acquire a T-dependent mass term are the 0o- They are called static modes 
because they are precisely those field components that do not fluctuate in T- 
direction. They are the only degrees of freedom that survive dynamically on 
scales L > T~^. All the other modes decouple. We call this scenario complete 
dimensional reduction. 

So far, the consideration is done on the level of the action, that is, for the 
classical field theory. How does this picture change under quantization? Now 
we have to investigate in full the partition function given by 



Z = 



Jv^c^p{-S{<P)) ■ Vc^^\{dc^{z), (58) 



with the action S{4)) as in (55). From the Matsubara-decomposition (56) it 
appears that the finite temperature field theory if the 0(2) in D dimensions 
can be viewed as a I? — 1-dimensional field theory of the ^„ (z) at zero temper- 
ature. There are multiple fields </)„ now almost all of which become massive. 
The only effect of the temperature is to renormalize the interaction. 

One might guess that on the infrared scale the heavy modes decouple as 
in the classical case, the only degrees of freedom left being the static fields 
(/)o(z) = (j)st{z). There are well known statements in quantum field theory on 
large mass behaviour. The Ambj0rn-Appelquist-Carazonne theorem predicts 
the decoupling of massive fields from the low-energy scale of the theory un- 
der quite general conditions. (It was proposed originally by Appelquist and 
Carazonne (1975). The proof of this theorem uses the BPHZ renormalization 
scheme and is due to Ambj0rn (1979).) In many cases it applies, such as for 
the simple scalar field theory above. The statement is very specific for the 
renormalization scheme used. All the coupling constants of the model have 
to be defined on the low momentum scale. In other instances the theorem 
does not apply. This concerns models with an anomaly or those subject to a 
non-linear symmetry, in which case the correlation functions are constrained 
by Ward identities. 
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The obstruction that prevents the appUcation of this theorem in our case 
is pecuhar. The number of heavy modes is infinite, whereas the above theorem 
does not make any prediction uniform on the number of modes. Anyway it 
would be a surprise if this theorem apphes to field theories at finite tempera- 
ture. For there is a strong interplay between renormalization and decoupling 
properties. Its validity would allow to circumvent the issue of ultraviolet 
renormalization of the theory, which is an intrinsic D-dimensional problem. 

Whether or not the quantum fluctuations destroy complete dimensional 
reduction must be decided by other means. The answer was given by Lands- 
man (1989) by means of a careful perturbative renormalization group anal- 
ysis. Apart from some exceptional cases (such as QED), the answer is nega- 
tive. To understand this fact, we have to account for renormalization theory. 
Let i^(p, T) denote a generic connected or 1-particle irreducible correlation 
function, p a generic spatial momentum. F achieves contributions both from 
the static and the non-static modes. The analogy of the classical decoupling 
property of the latter reads 

rip,T) ^ nt(p,T)(l + 0(|,^)), (59) 

where Fst denotes the purely static contribution to F. fi denotes some finite 
mass parameter. For renormalizable field theories we can always achieve this 
in the first instance, by imposing normalization conditions on some low mo- 
mentum scale. These conditions define the physical, renormalized coupling 
constants of the theory. The bare parameters of the model, that is the cou- 
pling constants of the action, are functions of the physical coupling constants 
and the ultraviolet cutoff. They are defined in such a way that the correlation 
functions F stay finite in the large cutoff limit and become functions of the 
renormalized coupling constants and the temperature only. This is the issue 
of renormalization. 

The point to be made now is that the bare parameters do not depend on 
the temperature T. The latter only enters through the torus of length T^^. In 
order to obtain (59), the subtractions must be temperature dependent, and 
so are the renormalized parameters. In particular, the renormalized mass 
/x that enters (59) becomes of the order of T in the generic case. Hence 
the correction term becomes of the order of 1. This shows that complete 
dimensional reduction breaks down. 

In accordance to the decomposition (56) of the fields into static and non- 
static components, let us write the partition function as 

V(f>st CXp{-Seff{c^st)) (60) 

with 

exp(-5e//(<^st)) = fV(bns exp{-S{T'^^(^st+(bns)) (61) 
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obtained by integrating out the non-static modes. Formally we can write this 
as 

Seff{(t>st) = S(3{4)st) + {6nsS){(t>st), (62) 

where Sq is the classically reduced action of the temporal zero modes (pat- 
This action is completely super-renormalizable. For the example above, it is 
given by 

^o(<^.*) - /^''"'^ (^cj}st{^){-AD-i4>st){^) + V{(l)st{^))\. (63) 

The quantum fluctuations of the non-static modes are completely described 
by the correction term 5nsS. This induced interaction is non-local. We know 
from the discussion above that it cannot be neglected even at very high 
temperature. Nevertheless one would still think of an appropriate high tem- 
perature expansion. Not any expansion will work because we have to avoid 
the same obstruction leading to (59). 

The solution is to generate the high temperature reduction as the low 
momentum expansion (Reisz (1992)). S^sS is made up of purely non-static 
modes. Thus it has an intrinsic mass gap. Its momentum representation is 
analytic at zero momentum p. We are thus allowed to expand about p = 
0. In a renormalizable field theory, the high temperature behaviour of the 
renormalized n-point function /«" is proven to be 

Tt\v,T) = TP-P,J^) + o(l) as T ^ ^. (64) 

Pm denotes a (multidimensional) polynomial of degree m for tti > 0, Pm = 
for TO < 0, and 

p„ = {D-l)-n^^. (65) 

The degree p„ depends only on n (and _D), but is independent of the order 
to which -Tn" has been computed in a perturbative (e.g. weak coupling) 
expansion. 

Let us write Tp/(p) for the (multidimensional) Taylor expansion about 
p = to the order 5 of the function /(p). To keep the finite T contributions 
on the infrared scale amounts to replace for those n with Pn > 

f(:)(p,T) ^ T^--''-r!,:\p,T) (66) 

and put all others to 0. Here, g„ denotes some integer number with < g„ < 
Pn- We notice that the expansion does not generate ratios of the form p-fT. 
The expansion is applicable even for T-dependent normalization conditions 
where the renormalized masses p increase with T. 

Only non-static correlations Pns with Pn > contribute to 5ns S at high 
temperature. The contribution is polynomial in momentum space, that is, 
Pns generates a local n-point interaction of SnsS. Hence we end up with a 
locally interacting model of the fields </ist(z)- 
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If we put Qn — above, i.e. keep the complete high T part of the non-static 
amphtudes, the effective model becomes a (power counting-) renormalizable 
field theory in D — 1 dimensions. 

We may be interested to keep only the leading term for every correlation 
function. In this case we put (/„ = /?„ in (66) and obtain 

fi:) (p, T) ^ T^r,t^ (P, T) ^ f (;') (0, T) (67) 

for all amplitudes Fn with /?„ > 0. This restriction is meaningful at suffi- 
ciently large spatial separations compared to the inverse temperature T~^, 
for instance for the investigation of universal critical behaviour at a second 
order phase transition. Corrections are suppressed by the order of p/T. 

We should emphasize at this point that (67) is not equivalent to keeping 
only those effective interactions that stay supcr-renormalizablc (except for 
the reduction from D = 3toD — 1 = 2 dimensions, cf. below). In general, 
in addition this would require to omit the contributions from the non-static 
correlation functions Fns with /?„ = 0. For example, for the reduction from 
4 to 3 dimensions this would imply to ignore interactions of the form 0®. 
This is justified only if the temperature is sufficiently large. Without further 
justification this omission is a potential danger for the investigation of the 
finite temperature phase transition. 

If only the effective interactions mediated by the non-static modes are kept 
which do preserve super-renormalizability of the (D-I)-dimensional effective 
model, the latter is made ultraviolet finite by the D-dimensional counterterm 
of the original theory, of course projected in the same way as in (67). In 
the more general case, supplementary normalization conditions have to be 
imposed. They are obtained by matching conditions to the D-dimcnsional 
theory. 

Let us summarize the state of the art for a local renormalizable quantum 
field theory at finite temperature. Its infrared behaviour on spatial length 
scales L > T~^ is quantitatively described by a local, generically renormal- 
izable or even super-renormalizable field theory of the purely static modes 
(/)st(z). This is a field theory in one less dimension and at zero temperature. 
The non-static degrees of freedom (j)ns, which account for the oscillations 
along the temperature torus, have an intrinsic mass gap of the order of T and 
hence do not survive as dynamical degrees of freedom in this region. However, 
they induce local interactions of the static fields that inevitably have to be 
taken into account in order to describe the infrared properties in a correct 
way. This concerns both the properties of the high temperature phase as well 
as finite temperature phase transitions. From the point of view of statistical 
mechanics, precisely these interactions are responsible for the phenomenon 
that quantum critical behaviour can be different from the critical behaviour 
of the corresponding classical model. Which interactions actually have to be 
kept depends on the situation at hand, as explained above. 
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A remarkable property of dimensional reduction is the fact that fermions 
do not survive the reduction process as dynamical degrees of freedom. As dis- 
cussed in the last chapter, the anti-periodic boundary conditions imposed on 
the fermion fields in the T-direction provide a mass gap, identifying fermions 
as purely non-static modes. The "only" effect they have on the infrared scale 
is mediated by renormalization of the interaction of the otherwise bosonic 
fields. 

The actual computation of the effective action of the static fields (pst can 
be done by various means. For four-dimensional gauge theories, QCD and the 
electroweak standard model weak coupling renormalized perturbation theory 
has been used extensively. For the Gross-Neveu model below we will carry 
out the dimensional reduction by means of the large N expansion. 

3.2 Application to the three-dimensional Gross-Neveu model at 
finite temperature 

We are back now to the Gross-Neveu model. This model is renormalizable 
for a large number N of flavours. Hence, dimensional reduction as presented 
in the last section applies. We have D = 3, so by (65) we get /9„=2 for all 
n. All non- vanishing correlation functions Fns {p,T) grow at most as T^ at 
high temperature. 

We are interested in the parity restoring phase transition and apply the 
high temperature expansion in the form (67) to all non-static effective ver- 
tices. That is, we keep the leading high temperature behaviour of /«" (p, T) 
for all n. It is convenient to rescale the twodimensional fields by a factor of 
(167r)^/^//3. Then the effective action is given by 

&„(«=/|^(i"*w^+^i:(^# 

rn>2 '' ^ 

Here and in the following we write instead of (\)st for the static fields. Only 
even powers of occur in the effective action. The bare mass squared TOq and 
the constants Cm for m > 2 are given by 

, 161n2 ,^ ^, to2 ^, 1 , 

(/3c - /3) + ^3277 + 0{ — ) 



/32/3^ '^'^ ' ' /32iv '7V2' 



■^..-nr^fe^ae'ra + oi-))- m 



The Cm do not depend on the number of flavours. Also, they are temperature 
independent non-negative numbers because of < Jm(0,/3) = const Z?^™""^, 
with ?7i-dependent constant. In particular, we have 
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C2 = 


= Stt, 




C3 = 


32 2 




C4 = 


128 3 
75 



(70) 



If A'^ = 00, only the quadratic part of the action survives and the model 
is purely Gaussian. On the other hand, for large but finite N, which we are 
interested in, we obtain an interacting model. The interactions are alternating 
in sign. For every m > 2, the 2TO-point coupling constant is of the order of 
jYi-m/2_ Pq], large N interactions (/)(z)^™ with large m are highly suppressed. 
The Z2-symmetric model described by the action (68) is super-renormaliz- 
able. The 1/N expansion is an expansion in the number of loops. The only 
ultraviolet divergence occurs for the 2-point function to the one loop order, 
which is the order A'^"^. There is only one corresponding, logarithmically 
divergent, mass counterterm included. This is the second term for TOq of 
(69). Its precise value does not play a role for the universality class of the 
phase transition, as will be seen in the next chapter. All other correlation 
functions stay ultraviolet finite and arc not subject to subtractions. 

4 Phase transition for finite number of flavours N 

We are prepared now to investigate the parity restoring phase transition of 
the Gross- Neveu model. This will be done by determining the phase structure 
of the associated 2-dimensional effective model described by the action (68) . 

Two-dimensional field theories typically evolve a rich and complicated 
phase structure. For instance, renormalization group studies of hierarchical 
models reveal that Z2 symmetric scalar models have complicated fixed point 
manifolds in field space (Finn, Pordt and Wieczerkowski (1994)). We have 
to expect more than just one (Ising) universality class for 2d scalar models, 
depending to a large extent on their interactions and the coupling constant 
strengths. 

The interactions and the values of the coupling constants of the effective 
model (68) have been computed by dimensional reduction. They are given 
functions of the number N of flavours and of the temperature T = (3~^ . The 
model is superrenormalizable, and there is only one logarithmically cutoff 
dependent mass counterterm. This is the only term that depends on the 
particular cutoff chosen. For the following non-perturbative investigation we 
put the model on the lattice. 

First of all, we truncate the action beyond the interaction (f)^. Terms of 
the order of N~'^ are omitted. This is legitimate for finite, but large N, which 
we are interested in. On the lattice, the action becomes 




xe/i2 ^ \ M=i,2 



I J2 (^M'/')(x) 
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Here /I2 denotes the two-dimensional hypercubic lattice with lattice spacing 
a, and d is the lattice derivative 

(9^0) (x) == - {cj^ix + aTi) - H^)) , (72) 

a 

fi as unit vector in the positive /i-direction. 

The choice of the lattice cutoff is consistent. For, the lattice model is 
(super-)renormalizable. According to well known power counting theorems 
on the lattice (Reisz (1988)), the continuum limit where the lattice spacing is 
sent to zero exists to all orders of the 1/N expansion. That is, all correlation 
functions of the basic field (j) stay ultraviolet finite. Furthermore, in this large 
cutoff limit they agree with those of the model with e.g. the simple momentum 
cutoff, which we have chosen so far. 

For further use we write the action in a form that is well known for lattice 
spin models. Furthermore, the lattice spacing is put to 1 in the following. 
We rescale the fields by introducing the so-called hopping parameter k and 
define further ultralocal coupling constants A, a and lu by 

(/)(x) = (2«)i/2 0o(x), 
k(4 + m^) = 1 - 2A + 3cr - 4a;, 

y2 

4c2K^— = A-3cr + 6u;, (73) 



,2-_ 
,T2 



-8c3K^— = a -Au, 

y2 
WcaK T — LO. 

Then, up to an irrelevant constant normalization factor, we obtain the par- 
tition function in the form 



Z = / n exp(-5o(0o)) (74) 

with the action 

So{M)= J2 ( 5^"H'/'o(x)) - (2«) J2 M^)M^ + n)], 

^(")(0o(x)) = 0o(x)2 + A(0o(x)2 - 1)2 (75) 

+a(0o(x)2-l)'Vc^(0o(x)2-l)'\ 

The coupling constants A, a and uj are of the order of A^^^, N^'^ and N~^, 
respectively, and are given by 
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A^36(l ' ''' 



with 



7^ 257^V ' 
1728 1 - 



= r..r and 7^ = - — . (77) 

9iV 27rK ^ ' 



The particular values where simultaneoulsy / = and TZ = oo correspond 
to the case of infinite number of flavours, TV = oo. The model then is a free 
Gaussian model. For all other values we introduce the ratios 

^ A 25^2 f 4 288 \ 

"^^ = 48^ V-n + 2^)^ 

^ o- 25 / 144\ 

As we will see, the phase structure can be parametrized as function of the 
ratios a and (3. 

The finite temperature phase transition and its properties will be deter- 
mined non-perturbatively by means of convergent high temperature series. 
We will apply the linked cluster expansion (LCE). It provides hopping pa- 
rameter series for the free energy density and the connected correlation func- 
tion of the model about completely disordered lattice systems. It should be 
emphasized that LCE generates convergent power series, convergence being 
uniform in volume. This implies that correlations can be measured with ar- 
bitrary precision within the domain of convergence, the precision depending 
on the order of computation. Even more, under quite general conditions, 
the phase transition is identified with the radius of convergence. The critical 
behaviour of the model is encoded in the high order coefficients of the se- 
ries. This concerns both the location of the transition as well as the critical 
exponents and amplitudes. 

The hopping parameter series have been generated recently for two-, four, 
and six-point functions up to the 20th order, for arbitrary coupling constants 
(Reisz (1995a)). The series are available for general 0(N) scalar field models 
on hypercubic lattices in £> > 2 dimensions. The lattices are unbounded in all 
directions or they have one compact dimension, representing the temperature 
torus. Also, the expansion has been generalized to finite lattices in order to 
allow for the determination of the order of a transition by finite size scaling 
analysis (Meyer-Ortmanns and Reisz (1997)). This LCE technique allows for 
a rapid generation of the series representation for all coupling constants. We 
will outline the idea behind this in the next section. 
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4.1 Convergent series expansion - LCE 

Hopping parameter expansions of classical lattice spin systems have a long 
tradition in statistical physics. Good references are Wortis (1974) and Itzyk- 
son and Drouffc (1989). In contrast to generic perturbation theory about 
eigenstates of harmonic oscillators or free fields, which arc normally at best 
asymptotically convergent, such expansions provide convergent series about 
disordered lattice systems. 

Linked cluster expansion (LCE) is a hopping parameter expansion ap- 
plied to the free energy density, connected correlation functions and suscep- 
tibilities. Providing absolutely convergent power series, under quite general 
conditions they allow for the determination of the critical properties of a sta- 
tistical model. This includes critical couplings as well as critical exponents 
and amplitudes. These numbers are encoded in the high order behaviour of 
the series coefficients. Their determination requires a computation to high 
orders. 

The realization of such expansions by convenient algorithms and with the 
aid of computers have been poineered by Liischer and Weisz (1988). They 
were the first to perforin such expansions to the 14th order of the hoppping 
parameter. Recent progress and generalizations have been made in various 
directions. They have been generalized to field theories at finite tempera- 
ture respectively quantum spin systems, where one dimension is compact 
with torus length given by the inverse temperature (Reisz (1995)). To the 
finite volume, allowing for finite size scaling analysis in order to determine 
the nature of a phase transition (Meyer-Ortmanns and Reisz (1997)). In the 
course of this construction it was possible to prolong the length of the scries 
to the order 20. Other progresses concern classical fixed-length spin systems 
(Butera and Comi (1997)) or higher moments of two-dimensional scalar mod- 
els (Campostrini, Pelissetto, Rossi and Vicari (1996)). 

Below we only give a short outline of the ideas behind the expansion 
technique. Details will be given only to the extent they are needed for the 
application in the next section. 

Let A]j denote a D-dimcnsional hypercubic lattice. For simplicity we work 
with units where the lattice spacing becomes 1. We consider a class of 0(N) 
symmetric scalar fields living on the sites of the lattice. The partition function 
is given by 

Z{H,v)^ I Y[ d^^{x) exp(-5(<Z',t;)-h^i/(x)-<Z'(x)). (79) 

xeAo xeA 

Here <1> — ('Pi, . . . ,'Pn) denotes a real, TV-component scalar field, H are ex- 
ternal (classical) fields, and 

N 



H{x)-P{x) = ^i7„(a;)<l'„(x). 



a=l 



N 
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The action is given by 

Sl-P,!;)^ ^ 5°(<?(x))-i Y. E <^a(xKfc(x,y)<Z',(y). (80) 

The "ultra-local" (0-dimensional) part S*" of the action describes the sin- 
gle spin measure. We assume that it is 0{N) invariant and that it ensures 
stability of the partition function (79). For instance, 

S"{$) = <P^ + \{<P^ - 1)2 + cr(<Z>2 _ 1)3 + ^(^2 _ j^^4^ (-g-^^ 

with Lu > 0. Dependence on the bare ultra-local coupling constants A, a and 
Lo is suppressed in the following. Fields at different lattice sites are coupled 
by the hopping term v{x, y). For the case of nearest neighbour interactions, 

, X r f 2k , X, y nearest neighbour, ,„_, 

^'^^(^'^^ - '^-^ I ; oirwise. ^^2) 

5 is the Kronecker symbol. 

The generating functional of connected correlation functions (free energy) 
is given by 

W{H,v)^\nZ{H,v), (83) 



Q2n 



dHa^ixi) ■ ■■dHa2„{x2n) 



W{H, v) 



■ (84) 

H=0 



In particular, the connected 2-point susceptibility X2 and the moment /i2 are 
defined according to 

X 

Kb /^2 = E ( E ^n < ^a{x)MO) >' . (85) 



X \ i—O 



The Fourier transform W{p) of the 2-point function is defined by 

^ (2^ e^''-(---) W{p). (86) 

The standard definitions of the renormalized mass mji (as inverse correlation 
length) and the wave function rcnormalization constant Zfi are by the small 
momentum behaviour of W, 

W{p)-' = ^iml+p^ + 0{p^)) asp^O. (87) 

Zr 
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From (85) wc obtain the relations 

2 

ml = 2D— , Zr = 2D^. (88) 

M2 M2 

The critical exponents ^,v^ri are defined by the leading singular behaviour 
at the critical point Kc, 

lnx2 =^ -7ln(Kc - k), 
lnm|j ~ 2v\n{Kc — n) , as k / Kc, (89) 

In Zji ~ z/T] In (kc ^ /i) 7 

such that lyrj = 21/ — 7. 

Let us assume for the moment that the interaction part (82) of the action 
between fields at different lattice sites is switched off, i.e. v = 0. Then the 
action decomposes into a sum of the single spin actions over the lattice sites, 
S{'P,v = 0) = X^a; '^°(^(^))' ^^^ t^*^ partition function factorizes according 
to 

Z{H,v = 0)= Yl Z^{H{x)). (90) 

xeAo 

It follows that 

W{H,v = Q)= Y, W\H{x)), (91) 

xeAo 

with W^ = InZ'^. In particular, the connected correlation functions (84) 
vanish except for the case that all lattice sites agree, i.e. for xi = X2 = ■ ■ ■ = 
X2n- In that case, 

W^i'l.„(2;i,...,X2„) = (2j!"i)!! ^2n(ai,...,a2„). (92) 

The C2n are the totally symmetric 0(N) invariant tensors with C2n(l, . • • , 1) 
= (2n — 1)!! = 1 ■ 3 ■•■... ■ (2n — 1). We have introduced the vertex couphngs 
as defined by 

^-- MF ^°(^ 



(93) 

H=0 



The vertex couplings W2„ depend only on the ultralocal part 5''^ of the ac- 
tion, that is on the single spin measure. They are obtained from the relation 
W^(H) = In Z"{H) or 

n>l ^ '' n>l ^ '' 

where 

o _ /d^J>J>fexp(-^"(g>)) 

'"" /d^<Pexp(-50(^)) ' ^^^^ 
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by comparing coefficients of the Taylor expansions about H = 0. The result- 

o^ o 

ing relations of the {i'2n} ^^'^ the {w2n} is one-to-one. To some extent, the 

o 

f 2n encode universality classes, as we will see below. 

The linked cluster expansion for W and connected correlations VK^^") is 
the Taylor expansion with respect to the hopping couplings v{x, y) about this 
completely decoupled case. 



W{H,v) = I exp^^«,b(x,y)— 4-— I W{H,v] 
1 ^^ dvab{x,y) I 



• (96) 



The corresponding expansions of correlation functions are obtained from (83) 
and (84). Multiple derivatives of W with respect to v{x,y) arc computed by 
the generating equation 

d ^^^,^^, 1 / d'^W dW dW \ 

-W{H) = - ,^,,,,^,,.,, + ^^,^^^^,...^ • (97) 



dvab{x,y) ' ' 2\dHa{x)dHt,{y) dHa{x) dHtiy) 

The expansion (96) is convergent under the condition that the interaction 
v{x,y) is sufficiently local and weak (cf. Pordt (1996) and Pordt and Reisz 
(1996)). Convergence is uniform in the volume. In all that follows we restrict 
attention to the case of (82), that is to euclidean symmetric nearest neighbour 
interaction of strength k. By inspection, for the susceptibilities such as (85) 
the coefficients of the power series of k are of equal sign. This remarkable 
property implies that the singularity closest to the origin lies on the positive 
real K-axis. This singularity then is the phase boundary or the critical point. 
In turn it is possible to extract quantitative information on behaviour even 
close to criticality from the high order coefficients of the series. 

We know that e.g. the 2-point susceptibility (85) is represented by the 
convergent series expansion 

X2 ^ Y. ^^^'2 (2«)^ (98) 

L>0 

with real coefficents a^^- From the high order behaviour of the coefficients of 
the series we extract the critical point Kc and the singular critical behaviour 
of X2 close to Kc- As 

X2 ^ ("l - ^"j , (99) 

we have 

^^ = ^(l + ^+o(i^-)). (100) 

aL_i,2 2kc \ L J 

Analogous relations hold for the other correlations, as for the correlation 
length m]^ and for Z^ . Similarly, critical amplitudes are extracted. There 
are more sophisticated and also more involved methods than (100) to obtain 
the critical numbers to high precision. Our application below is to distinguish 
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Ising and mean field behaviour in two dimensions, which have rather different 
critical exponents. For this the ratio criterion (100) turns out to be sufficiently 
precise. 

The explicit computation of the series becomes rapidly complex with in- 
creasing order. It is mainly a combinatorical problem and is conveniently 
handled by graph theoretical terchniques. Correlation functions are repre- 
sented as a sum over equivalence classes of connected graphs, each class being 
endowed with the appropriate weight it contributes. For instance, the coef- 
ficients aL.2 are obtained as the sum over the set Gl,2 of connected graphs 
with 2 external and L internal lines, 

regL,2 

with w{r) the weight of the graph F. Let Vr denote the set of vertices of F, 
and for every v e Vr let l{v) be the sum of internal and external lines of F 
attached to the vertex v. The weight u'(r') then is given by 

^.(^) = 7^(n.^ (5|f^- (102) 

The first factor Tl{F) is a rational number. It is the product of the (inverse) 
topological symmetry number of F, the internal 0(N) symmetry factor, and 
the lattice embedding number of the graph F. This is the number of ways 
F can be put onto the lattice under the constraints imposed by the graph 
topology, appropriately supplemented with possible geometrical weights, as 
e.g. for /i2 of (85). Two vertices have to be put onto nearest neighbour lattice 
sites if they have at least one internal line in common. Beyond this there are 
no further exclusion constraints. In particular, a lattice site may be convered 
by more than just one vertex (free embedding). 

The second factor of (102) is the product over all vertices u of _r of the 

vertex couplings Vuy) (5**^). This is the only place where the ultralocal cou- 
pling constants of S*" enter. Generally this factor is a real number, it becomes 
rational only in exceptional cases. An example of a weight w{F) is given in 
Fig. 5. 




{^)i2Dfi9{N + 2)iN + 4)) (iS:)'(i^:) 



Fig. 5. A graph contributing to the connected 2-point function of the D-dimensional 
0(N) model, and its weight. 
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For given dimension, lattice topology and internal symmetry, TZ{r) is a 
fixed rational number. Its computation is expensive, but needs to be done 
only once. Changing the parameters of the action (that is, k and S*") does 

not require to compute TZ{r) a new. S^ only enters the vertex couplings W2„. 
Once we have computed all the Tl{r), the computation of the LCE series for 
a new S*" requires only a minimal amount of computing time. 

The number of graphs F and in turn of the Tl{r) is large. To the 20th 
order there are more than 10^ numbers TZ{r). Actually, it is not necessary 
to keep the T^{r) for all graphs F of interest. The idea to circumvent this 
is to introduce the so-called vertex structures. Roughly speaking, a vertex 
structure accounts for the numbers of vertices with particular numbers of 
lines attached. More precisely, a vertex structure w is a sequence of non- 
negative integers, 

w = {wi,W2,...) , w.G {0,1,2,...}, (103) 

with only finitely many Wi ^ 0. We write W for the set of all vertex structures. 
For a graph F the vertex structure w = w{F) G W associated to F is defined 
by (103) with Wi equal to the numbers of vertices of F that have precisely i 
lines attached. With the definition 

we obtain instead of (101) and (102) 

(oc \ -Wi 
(7^ • ^'^'^ 

The sum is finite because most of the TZ{w) vanish. Actually there are only 
a couple of hundreds of non-trivial 7?.(w)'s, to be compared to 10^ TZ{Fys. 

In order to make high orders in the expansion feasible it is necessary to 
introduce more restricted subclasses of graphs such as 1-line and 1-vertex ir- 
reducible graphs and renormalized moments. The correlations are then repre- 
sented in terms of the latter. The major problems to be solved are algorithmic 
in nature. Examples are: 

— unique algebraic representation of graphs 

— construction and counting of graph 

— computing weights, in particular lattice embedding numbers. 

High orders are feasible by largely exploring simplifications. On the hypercu- 
bic lattice, only graphs matter that have an even number of lines in each loop. 
For 0(N) symmetric models, vertices with an odd number of lines attached 
vanish. The most important series are available to the 20th order yet. 
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4.2 The nature of the parity transition 

We determine the phase transition of the model (74), (75). Before we put the 
coupling constants to their prescibed values given by (76), it is instructive to 
discuss the strong coupling limit first. 

As we have mentioned in the last section, universality classes can be read 

o 

off to some extent from the n-dependence of the vertex couplings W2n as 
defined in (95), without having to go through the complete analysis of the 
LCE series. For Z2-symmetric models, the two cases we are interested in are 
the following. Let z be any positive real number. Then 

s,,.ei_i^,. (100) 

for all fc = 0, 1, 2, . . . implies that the model belongs to the universality class 
of the Gaussian model. The critical hopping parameter Kc is given by k^z — 
1/{2D) in D dimensions, and the critical exponents are 7 = 1, i^ = 1/2 and 
?7 = 0. On the other hand, let the vertex couplings be given by 

o _ (2fc-l)!! r(i) 

Then the universality class is that of the Ising model. In two dimensions, 
KcZ = (1/4) In (2^/2 + 1), 7 = 1.75, u = 1 and t] = 0.25. 

The strong coupling limit of (74), (75) is defined by w ^ 00 with the 
ratios a ~ X/lo and /3 = a/uj hold fixed. The behaviour of the vertices V2k is 
obtained by the saddle point expansion. As function of a and P the model 
evolves a rather complicated phase structure. We don't need to discuss it in 
detail, but the following regions are of particular interest to us. 

There are regions where the behaviour is Gaussian, and regions where it 
is Ising like. For S//3^ > 1/4 we obtain (107) with 2; = 1. In this case, Kc = 
(1/4) In (2^/2 + 1) independent of a and /3. For < a//?^ < 1/4 and /3 < we 

again get (107), with z = 1- (3/3/8) ( 1 + (l - f ^) ) > 0. In both cases 

the transition is Ising like. On the other hand, in the region < a//3^ < 1/4 
and /? > 0, S < min (|/3 — 2, /? — 1) implies (106), with z ^ oj^^. The model is 
Gaussian behaved or completely decouples over the lattice because of z ^ 0. 
In Fig. 5 we show a qualitative plot of the universality domains. 

In our case the ratios a and (3 are largely fixed. They are given by (78), 
and 

^ ov 1 4 I 288 

P^ 25 (i_^)^ 

Because TZ is proportional to N, for sufficiently large N we have both /3 < 
and a//3^ > 1/4. This inevitably selects the Ising universality class. 
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Fig. 6. Part of the phase structure of the effective model in the strong coupling 
limit. In the regions I the critical exponents are those of the 2-dimensional Ising 
model. Region II is the Gaussian domain or decoupled across the lattice. 



Next we consider the case of finite coupling constants A, a and w. We 
invoke the LCE to construct the large order series of the susceptibilities X2 
and ii2 as defined in (85) for the various coupling constants. From the series 
coefficients then we extract the critical coupling Kc and the critical exponents 
V and 7, according to (99) and (100) and in analogy for the correlation length 

We proceed in the following way. We choose TZ — 0{N) sufficiently large. 
Then a and j3 are fixed. For every TZ we vary the 8-point coupling lo from ui = 
oo down to small values, the other couplings running according to A = aui and 
a = [3uj. From the LCE series of X2 we obtain the critical hopping parameter 
Kc(A, CT, w). This then determines the precise value of N corresponding to the 
transition point by 



2ttKc ^ 

N = — - n. 



(109) 



Then the critical exponents 7 and v are measured. As a result, their values are 
independent of w, and they are equal to the Ising model ones. The exception 
is a small region about the Gaussian model with A = cr = w = 0, where there 
is a smooth crossover to the Gaussian numbers. This is to be ascribed as a 
truncation effect of the finite order series. The above process is repeated for 
various values of TZ so that N varies over a large region. 

In Table 1 and Fig. 7 we have collected some numbers on the critical 
exponents 7 and v obtained in this way. For all N the critical exponents 
are those of the 2-dimensional Ising model. For small N this might be not 
trustworthy because we have obtained the effective model within the large N 
expansion. On the other hand, for large N this is predictive. 

As a result, we conclude that the 3-dimensional Gross-Ncveu model, with 
large but finite number A'^of fermionic flavours, has a finite temperature phase 
transition that belongs to the universality class of the 2-dimensional Ising 
model. 
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N 


7(iV) 


2iy{N) 


46.1 
24.3 
4.68 


1.751(2) 
1.752(3) 
1.749(2) 


2.001(2) 
2.001(4) 
1.999(2) 


1.41 


1.750(2) 


2.001(2) 



Table 1. The critical exponents 7 and v of the three-dimensional finite temperature 
Gross-Neveu model for various numbers A'' of flavours. They agree to those of the 
two-dimensional Ising model. 



5 Summary 



We have investigated the finite temperature phase structure of the three- 
dimensional Gross-Neveu model. It is a parity symmetric, four-fermion inter- 
acting model with A^ fcrmion species, and we kept N large but finite. 

From the point of view of high energy physics, the interest in this model 
is based on the observation that it reveals remarkable similarities to the 
chiral phase transition of QCD with two massless fiavours. Being a three- 
dimensional model, chirality is replaced by parity, but both symmetries would 
be broken explicitly by a fermionic mass term. 

Beyond its similarity to QCD, the advantage of this model is that to a 
large extent we have a closed analytic control of it. This includes the finite 
temperature transition. 

At zero temperature, parity symmetry is spontaneously broken. It be- 
comes restored at some critical temperature Tc by a second order transition. 
There are at least two candidates of universality classes this transition might 
belong to. These are the two-dimensional Ising model and mean field be- 
haviour. Both scenarios could be true, following different conventional wis- 
dom arguments. It became clear in this lecture that the definite answer can 
be given only by explicit computations. 

Our result is that for large but finite number TV of fermions, the transition 
belongs to the universality class of the two-dimensional Ising model, with 
critical exponents 7 — 1.75, v — I and 77 — 0.25. This result is obtained 
by the combination of three computational techniques: Large N expansion, 
dimensional reduction in quantum field theory and high order linked cluster 
expansions. 

The large N representation of the Gross-Neveu model amounts to a re- 
summation of planar graphs, that is, of chains of fermion bubbles, cf. Fig. 3. 
In the path integral approach this can be achieved in closed form by repre- 
senting the four-fermi interaction as the result of a Yukawa-type interaction 
with an auxiliary scalar field. Whereas the weak (quartic) coupling expansion 
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is not renormalizable, the model becomes (strictly) rcnormalizable for large 
TV, in particular in the expansion in powers of N^^. 

Dimensional reduction is to be considered as a special high temperature 
expansion that is applied in the course of integrating out the degrees of free- 
dom that fluctuate along the temperature torus. It is applied to a quantum 
field theory (rather than a classical field theory) to describe its infrared prop- 
erties. Applied to the Gross-Neveu model, the latter is mapped onto a two- 
dimensional effective field theory. Due to the anti-periodicity of the fermionic 
fields, the model is purely bosonic. It is local and super-renormalizable. The 
interactions of this model are computed within the large N expansion. Both 
because the transition at the critical temperature Tc is of second order and 
N is supposed to be sufficiently large, dimensional reduction will work down 
toTc. 

Finally, the relevant phase structure of the dimensionally reduced effective 
model is determined. This is done by applying the linked cluster expansion. 
This technique allows us to compute convergent hopping parameter series 
for the free energy and connected correlation functions to very high orders. 
It provides the series also for the correlation length and the wave function 
renormalization constant. The high order behaviour of the series coefficients 
allow for the determination of critical couplings, including critical exponents. 

The effective model has a complex phase structure. In particular there 
are both regions in coupling constant space with mean field and with Ising 
behaviour. However, the coupling constants are not free but have been com- 
puted already by dimensional reduction. They identify the Ising region as 
that part which describes the finite temperature Gross-Neveu model. 

It would be interesting to supplement these considerations by other means. 
In the large N expansion, N^^ acts as a real coupling constant. The parity 
symmetry of the model does not depend on N. This promises an investigation 
by means of the renormalization group, with N^^ kept as a running coupling 
constant. This investigation should reveal a non-trivial infrared fixed point 
at some finite N* < cx). For instance, one could try an e-expansion about 
1 + 1 dimensions. 

What do we learn for the QCD chiral phase transition? The situation 
is similar to the Gross-Neveu model. The effective model now is a three- 
dimensional, pure bosonic theory. Fermions do not act as dynamical degrees 
of freedom on the infrared scale R > T^^, but they determine the bosonic 
interactions. The cr-model scenario of the two flavour case is not ruled out, 
with critical exponents of the three-dimensional 0(4) model. This is also 
supported by recent numerical investigations (Boyd, Karsch, Laermann and 
Oevers (1996), cf. also the contribution of Frithjof Karsch to these proceed- 
ings). 

Dimensional reduction has been proved to be a powerful method to study 
infrared properties of the QCD plasma phase. The reduction step from four 
to three dimensions is normally done by means of renormalized perturba- 
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tion theory, that is as an expansion with respect to the rcnorniahzed gauge 
couphng constant g^-cniRjT) at temperature T and length scale R. In this 
way there is analytic control of the cutoff and volume dependence. On lat- 
tices with moderate extension in the temperature direction, these effects are 
rather large. This implies a considerable complication for numerical investi- 
gations. Their analytic control allows to remove these effects. In this way it 
is possible to determine screening masses of the quark gluon plasma already 
on small lattices. 

This works fine for temperatures at least above twice the phase tran- 
sition temperature Tc- Compared to the Gross-Neveu model we have one 
more spatial dimension. From dimensional power counting this implies that 
interactions of higher operators are suppressed by inverse powers of the tem- 
perature. We do not have to rely on a large number of flavours any more to 
truncate those terms. However, approaching the phase transition, the renor- 
malized gauge coupling constant becomes large, and we can no longer rely on 
perturbation theory. The relevant coupling constants of the effective model 
have to be computed by non-perturbative means. Therefore, the universality 
class of the two-flavour chiral transition in QCD is still an open question. 
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